Entanglement entropy of highly degenerate states and fractal dimensions 
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We consider the bi-partite entanglement entropy of ground states of extended quantum systems 
with a large degeneracy. Often, as when there is a spontaneously broken global Lie group symmetry, 
basis elements of the lowest-energy space form a natural geometrical structure. For instance, the 
spins of a spin-1/2 representation, pointing in various directions, form a sphere. We show that for 
subsystems with a large number m of local degrees of freedom, the entanglement entropy diverges as 
| logm, where d is the fractal dimension of the subset of basis elements with non-zero coefficients. 
We interpret this result by seeing d as the (not necessarily integer) number of zero-energy Goldstone 
bosons describing the ground state. We suggest that this result holds quite generally for largely 
degenerate ground states, with potential applications to spin glasses and quenched disorder. 
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The entanglement entropy is a measure of entangle- 
ment between two complementary sets of observables in 
a quantum system [IJ . It is denned as the von Neumann 
entropy of the reduced density matrix of the state \^>) 
with respect to a tensor factor of the Hilbert space H: 

S = -Tr A {p A \ogp A ) with PA = Tr B \V){V\, (1) 

and H = A <E) B. A related measure is obtained from the 
Renyi entropy, S n = lo § Tt APa)> clearly, S = S% = 
lim„_ ! . 1 + S n . The entanglement and Renyi entropies have 
important applications to e.g. quantum computation and 
numerical simulations of quantum systems. 

In extended quantum systems near to critical points, 
the entanglement entropy has turned out to reveal funda- 
mental properties of ground states (for reviews, see e.g. 
0). An important result is the so-called area law. Con- 
sider a quantum system of dimensionality D > 2 with 
correlation length £, and a subsystem A composed of the 
local degrees of freedom on a Z?-dimensional region A of 
linear extension I (generically, the region A is composed 
various components of different connectivities, and i is 
the overall scale of A) . It turns out that the entanglement 
entropy between the subsystem and the rest diverges as £ 
and £ increase, the ratio r — £/£, being fixed, with a power 
law £- D ~ 1 , withpossible logarithmic corrections for gap- 
less systems [3|-|8[ . But this area law is special in the cases 
where D = 1. There, the divergence is always logarith- 
mic: 2£ log(£) where q is the number of points separating 
A from the rest and c is the central charge of the crit- 
ical theory 0, [l(|. Interestingly, the number c comes 
out, which essentially measures the number of degrees 
of freedom that are carried over from the microscopic 
theory to the macroscopic universal theory. Further, for 
D = 1 again, subtracting this divergence, the rest is a 
finite quantity which depends on r, which saturates to a 
finite value at r — oo, and which tends to this value in 
an exponential way that is solely determined by the spec- 
trum of masses of the corresponding perturbation of the 



critical point llj, Il2j . The spectrum of asymptotic par- 
ticles characterizes the low-energy degrees of freedom of 
the universal theory. Moreover, in systems with a bound- 
ary, the boundary degeneracy appears also by a natural 
subtraction [H El. This de generacy characterizes the 
number of degrees of freedom carried by the boundary. 
These results point to the observation that if the entan- 
glement entropy diverges logarithmically at large subsys- 
tem size I, then the way it diverges is controlled by some 
basic counting of universal degrees of freedom. 

All these results were established for non-degenerate 
ground states, or ground states with small, finite degen- 
eracies. A question arises as to the entanglement entropy 
for highly-degenerate ground states. Let us start by dis- 
cussing an example where a symmetry group is sponta- 
neously broken: the Heisenberg ferromagnet. This is an 
iV-site lattice L with spin-1/2 local degrees of freedom, 
with Hamiltonian: H = J^fy) e edges ofh S i' S h J < 
(Si is a vector of Pauli matrices acting on site i). This 
model has an SU(2) global symmetry, and the states 
— <8>iei,|Vv)ij where all spins point in the same 
direction v (i.e. Si ■ v\tpv)i — \ipv)i, \v\ = 1), span the 
lowest-energy subspace. In the usual description, we 
make a choice of an arbitrary direction vq. Such a state 
is not invariant under SU(2) transformations, hence the 
symmetry is dynamically broken. The Hilbert space rlv 
in the thermodynamic limit N —} oo is then, with the 
ground state IVVo)^! the se t of all finite-energy, local 
excitations above it. By locality of the Hamiltonian, it 
excludes all ground states and excited states associated to 
other directions, rl$ for v ^ vq (these cannot be reached 
by a finite number of local changes of the infinite system) . 

But linear combinations of IV^^s also give lowest- 
energy states, and for them we will take a different de- 
scription of quantum states that is more appropriate. 
The lowest-energy subspace is the ./V+l-dimcnsional sub- 
space forming a spin-iV/2 representation. For every N, 
there exists a set of N + 1 points Vk such that the set 
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of vectors IV'ik)^ forms a basis for this subspace. Fur- 
ther, in the limit N — ¥ oo, every point on the unit sphere 
is arbitrarily close to such a basis point. A good de- 
scription of the resulting space of infinite- volume lowest- 
energy quantum states is then obtained by using the ge- 
ometry induced by averages of local operators (see [14J , 
Section 4) . In this geometry, the distance between states 
in directions v and v' is smoothly related to the distances 
between the vectors v and v' on the unit sphere. This 
geometry is convenient for discussing the entanglement 
entrop y, b ecause, as is developed in [14J (based on earlier 
works [ll|]), the latter can be evaluated from the average 
of a local permutation operator. 

Linear combinations could involve infinitely many di- 
rections v, with appropriate integration measures on the 
unit sphere. This occurs, e.g., when a ground state of the 
infinite-length one-dimensional Heisenberg ferromagnet 
is reached by an adiabatic lowering of the anisotropy of 
the XXZ model: an integration over a great circle on the 
unit sphere is obtained [14j]. Although each state l^v)^ 
has zero entanglement entropy (since it is factorisable) , 
linear combinations do not, and linear combinations in- 
volving infinitely many directions v should have growing 
entropy as I — > oo. What is the £ — > oo behaviour for 
such infinite linear combinations? 

Let A be composed of m degrees of freedom and 
N = oo. Clearly, any minimal-energy state has a sym- 
metry under exchange of any two sites, hence the entan- 
glement entropy depends on m but not on the particular 
sites chosen. We may take the m sites to form a contin- 
uum of dimension D, writing m = £ D . First note that a 
large-m divergence \ log m of the entanglement entropy 
was found in 14|, Il5f for the state formed by an integra- 
tion over a great circle. Second, recall that when there 
is spontaneous symmetry breaking, there are massless 
excitations in the spectrum, the Goldstone bosons. In 
general, our idea is that linear combinations composed 
of all points along an arc on the unit sphere should be 
interpreted as representing the presence of a zero-energy 
Goldstone boson corresponding to the continuous motion 
along this arc. Moreover, every linearly independent lo- 
cal direction on the unit sphere corresponds to a linearly 
independent Goldstone boson, each of which can be seen 
as a universal degree of freedom. Hence, the observations 
above suggest a divergence of the form | log m, where d 
is the number of Goldstone degrees of freedom present 
in the linear combinations. This number is simply the 
dimension of the support of the linear combination on 
the unit sphere. The "number" of Goldstone degrees of 
freedom d is not restricted to the integers: the support 
of the linear combination may have a fractal dimension. 

Here we argue that the large-m (large-€) behaviour is 



5„ = ^logm + 0(l) = ^log£ + 0(l) (2) 
for all n, where d is the (fractal) dimension of the support 



of the linear combination, < d < 2 for the spin-1/2 
Heisenberg ferromagnet. For instance, if the great circle 
in the above example is replaced by the Cantor set, we 
would find d = log(2)/log(3). Note that the result is 
independent of n. 

A simple application of this formula is to detect a pos- 
sible blurring of the dynamically chosen direction vq ob- 
tained, for instance, as the system is cooled in a fixed 
magnetic field. The blurring could lead to a linear com- 
bination covering a possibly fractal small surface around 
vq. Standard local observables would not discern this, 
whereas formula @ shows that the entanglement entropy 
is extremely sensitive to it, giving d > instead of d = 0. 

Beyond the Heisenberg ferromagnet, our derivation be- 
low makes it clear that © should hold much more gen- 
erally. There are two conditions: 1) there exists a basis 
for the lowest-energy subspace where the entanglement 
entropy of each basis element is zero or small, and 2) 
the basis elements are given the local-operator geometry 
14j . The dimension d in (2} is that of the support of 



the linear combination in this geometry, which may be 
an integer or not, and which essentially counts the num- 
ber of Goldsone bosons in the state. For instance, for a 
quantum system in a "Mexican hat" potential, the ba- 
sis set is the geometrical circle at the bottom of the hat, 
and subsets of this will have < d < 1; similar obser- 
vations hold for any system with spontaneously broken 
continuous symmetry. 

Note that some "permutation-symmetric" (PS) states 
of the above type in a Hilbert space with on-site spin S 
were considered in [la ] . and ([2]) with d = 2S was ob- 
served. An analysis shows that the states chosen gener- 
alize the 5=1/2 great-circle state. Further, this is in 
agreement with our general arguments, which imply that 
all possibilities < d < AS may occur for PS states. 

Our formula is in sharp contrast with behaviors re- 
viewed above (e.g. £ D_1 for D > 1), related to the ge- 
ometric structure of the region A and arising thanks to 
locality of the system's interaction. To explain this, con- 
sider the case of a spontaneous symmetry breaking: local 
interactions only fix the lowest-energy subspace, not the 
ground state. By choosing a basis of lowest-energy states 
with minimal entanglement and with a local-operator ge- 
ometry, we expect that we encode all local information 
in the basis states themselves, and only the symmetry 
is probed by the degeneracy. Hence, our result, which 
has to do with the degeneracy, cannot measure geomet- 
rical objects in the system's real space. Rather, a loga- 
rithm occurs, whose coefficient counts Goldstone degrees 
of freedom, associated with the symmetry. If condition 
1) above does not hold, we expect two contributions to 
the asymptotic of the entanglement entropy: that of the 
large degeneracy ([2]), and that coming from locality. 

Examples of fractal sets of minima are found wher- 
ever random potentials occur (quenched disorder, see 
e.g. [HI), with possible connections to glasses. In the 
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classical phenomenology 17|, at temperatures below the 
glass transition point, the free energy surface reveals finer 
structures in the form of new energy minima within previ- 
ous valleys, displaying self-similarity and a fractal struc- 
ture; the set of effective minima has a nontrivial fractal 
dimension. High classical degeneracy also naturally oc- 
curs in frustrated spin systems [l8j]. These degeneracies 
may be lifted by quantum fluctuations (so-called "order 
by disorder"), although the underlying classical degen- 
eracy is known to have nontrivial quantum effects and 
to survive semi-classically By our formula the 
entanglement entropy could provide a further indicator 
at the quantum level of this (semi-) classical degeneracy. 

In the rest of this paper, we provide the main lines of 
the derivation of (J2)). A more precise proof and statement 
will be given in a separate work. 



THE HEISENBERG FERROMAGNET CASE 

The present derivation uses the replica trick, whereby 
n is assumed to be an integer > 1. The result, how- 
ever, can be interpreted for n 6 (l,oo). This provides 
the unique analytic continuation which does not diverge 
exponentially at large n; as is usual, this analytic contin- 
uation is assumed to provide S n for all real n > 1. 

Given a point v on the unit sphere S 2 , let us de- 
note by ipv G J~ the quantum state corresponding to 
the N -> oo limit of |7/v) (Ar) - As developed in [14j . 
a quantum state is a linear functional on the space of 
finitely-supported operators, which evaluates the aver- 
age; e.g. ^{O) = lim^oo W<^|0|^>W. We can 
write ips as a product of single-site quantum states, all 
acting in the same way: ipjj = (^) igZ tpv-i ■ At infinite vol- 
ume, vectors pointing in different directions have zero 
overlap: Kmjv^oo {N) &v\ipj) {N) = for v ^ v'. This 
holds as well with insertions of finitely-supported opera- 
tors, so the infinite-volume limit of linear combinations 
^tfdvlipv)^ gives the quantum state 



:= E Kl 2 ^' E I" 5 ' 2 = L 



(3) 



In order to evaluate the Renyi entanglement entropy 
associated to the ground state ^{ a ^} we recall the ap- 
proach developed in [3]. There, the Renyi entropy of 
a region A in a quantum state ip was expressed as an 
average on the n th tensor power of i\>: 



S n — 



1 



■log ty® n (7A)) 



(4) 



The operator averaged is Ta — TiieA^,, where 71 are 
local cyclic replica permutation operators which act on 
site i of the quantum spin chain by cyclicly permuting 
the spins of the n replicas of the model at that particular 



site. One of the results of [14| was the closed formula 

n 

7- = Tr aux n Y, E ^ 2 KT, (5) 



a— 1 ei,62 



where Ey £l represent elementary 2x2 matrices with a 
single non- vanishing entry at row £2, column e\, acting 
on space V — a, i (site i tensor copy a) or V = aux 
(auxiliary space). For the quantum state ij)i a ,\, we write 



{a„-} 



e n 1- 

{v a :a=l, ...,«,} \a=l 




From the trace expression ([5]) we find 

<g) (Ta) = I] T W II E E l^v a ;i {E? C1 ) ■ 

a i£A a e±,e 2 

Clearly, ipv a -i (-Ej 2£l ) is independent of i. Writing \ipc) = 
Sv,\ \ t)+Su,2| I), wefind^.j (E^f 1 ) = s^ e2 s^ a<ei , and 
tracing over the auxiliary space we obtain 

Tr aux J] E££*l>*a;i ( E ?^) = l^- + i>- 

Oi 6l,62 

Hence, we find 



Sn = log f E 



n 



with v. 



n+l 



= v\ . This saturates at large m to 
lim S„ = _!— logf VK| 2 " ) . 



(7) 



(8) 



That is, as expected, for any ground state given by a 
finite linear combination of basic zero entropy states, the 
entanglement entropy reaches a finite maximum as the 
number m of site of A tends to infinity. This corresponds 
to the case d = in ([2]). 

More interesting behaviours are obtained from "infi- 
nite linear combinations" of basic states, generalising ([3]). 
Given a smooth, self-avoiding path 7 : [0, 1] — > S 2 and a 
smooth function / : S 2 -)• E+ with J* |df(t)| /(•?(*)) = 1, 
the following integral can be defined and is a quan- 
tum ground state: tp^ '■— Jq |^t(*)| f{l{t)) V^t*)- Sim- 
ilarly, let /2 : [0, 1] x [0, 1] — > S 2 be a two-dimensional 
smooth curve such that J J Q \d 2 fi(X, 0)|/(/2(A, <]))) — 
1 (where \d 2 jl(X, cf>)\ is the surface element on the 
unit sphere). Then, the following is a ground state: 

V" (2) : = Jo 1 Jo I^mO, 0)1 /(m(A, </>)) V^(A»- Generalising, 
we may consider the set of non-zero coefficients to be 
a fractal set W C S 2 with fractal dimension d. With 
dH(v) the corresponding Hausdorff integration measure, 
we may write 

^ {d) := [ dH{v)f{v)iJv with f dH{v)f{v) = l. 



hv 



Iw 
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The Hausdorff measure is expected to occur naturally in 
taking the large- volume limit, if the set of vectors v such 
that a# 7^ becomes a fractal set. 

Computing the Renyi entropy of ip^ yields a simple 
generalisation of (|7|) where the sums over the vectors v a 
are replaced by integrations and the coefficients \aa a \ 2 
by the functions f{v a ). The logarithmic factor in 
becomes 



log / H dH(v a ) f(v a ) JJ<^„ \^ a+1 ) 



w 



(9) 



The large to asymptotics of these expressions will how- 
ever be radically different from (jHJ: there will be no sat- 
uration, and we will recover the behaviour highlighted in 
(|5J). This can be shown using a saddle-point analysis, as 
was done in [lij ] in a particular case; here it is generalised 
to integrals over fractal domains. 

For the explicit calculations, we use = 

-j- ( ~ Z ie 1 ' wnere v =: (x, y, z) is a unit vector, 

and 9 is defined by x + iy = yl — z 1 e ie . From this 
we see that KV'tH^to)! < 1, with equality if and only if 
v = w. The saddle-point analysis from © is done by 
expanding the overlaps {ipv a \'4>v c , +1 ) around v a = v a +i, 
and re-writing the product Y[ a °^ these overlaps as an 
exponential. We get 



J[(' t PvJ^v a+1 ) = exp 



1 



n 

E 

a=l 



v a +i-v a \ + 



(10) 

where the ellipsis stand for terms that are order-2 and 
antisymmetric, and higher order terms. The order-2 an- 
tisymmetric terms vanish when the integrations in © 
are performed. 

In the case of tp^ 1 ' f° r instance, we may use the as- 
sumptions relating to / and 7: both are smooth, and 
the curve 7 is self-avoiding. Hence, with v a = 7(£ a ), the 
maximum occurs when t a — i a +i f° r all a = 1, . . . , 71. In 
this case we find 



exp 



a=l 



(11) 

where I7I 2 is evaluated at t = t\. The saddle-point 
analysis is then performed as follows. We need to raise 
the quantity above to the power m and substitute into 
the integral (j9]). We can replace f{^{t a )) by /(7*(ti)) 
for all a, since / is smooth. Changing variables to 
ti = \/m(ti — t\), i = 2,...,n guarantees that larger 
positive powers of t a give lower-order contributions at 
large m. We obtain 



S n — 



1 



log 



1 



dhmhTfm^y 



(12) 



He 



e (t a+ i-t a y+4+ti 

a = 2 



+o(t 3 /V^) 



where the integrals over ti , . . . , t n have been extended 
to (—00, 00) (the resulting correction terms are exponen- 
tially small). These integrals are of standard gaussian 
type and can be carried out explicitly. 

A very similar computation can be carried out for the 
state ip 1 " 2 ^ instead of ip^ 1 '. The final result can be ex- 
pressed in both cases d = 1 and d — 2 as 



d fm 
i„ ~ — log — 
2 6 V87r 



1-71 



loe 



\d d v\ f(vY 



(13) 

where higher order corrections would be O (to 1 / 2 ). This 
is in agreement with ([2]). 

Note that the constant term in (fTS")) is — log(/ max ) as 
n — > 00, where / max is the maximum of / on its support. 
At n — 1, we have rather d/2 — J \d d v\ f(v) log f(v). 

The calculation for fractal sets follows similar lines. 
A crucial feature of the Hausdorff measure is its scaling 
property. On the plane, the Hausdorff measure H' satis- 
fies s d H'(W) = H'(sW + u) for any W C R 2 . For the 
measure % on S 2 , this scaling covariance is replaced by 
an asymptotic behaviour that gives rise to the measure 
H' on the tangent plane: 

lim TO d/2 dH{il t /Vm + n) = dH'itii). (14) 

m—>oo 

Putting (fT0|) inside ([9]), changing variables to u; = 
V / ro(t/i — vi), i = 2, . . . , ti, and using ([M]). as to — > 00, 



1-71 



log 



d(n-l) 

TO 2 JW 



~ n 

dn(vt)f(vi) / [] 



<m'{vi)f I -j= + vi 1 e 



E l^ + l-^c,| 2 + K2| 2 + |ff„| 2 



where W m = v / to(W — v\) (projected to the tangent 
plane at vi). Although the integral might not exist in the 
large-TO limit, it is bounded, thanks to the exponentially 
decaying factor. This boundedness immediately gives rise 
to the leading asymptotics ©. It would be desirable 
to investigate more precisely the nature of the constant 
corrections to this general leading behaviour; we hope to 
return to this in a future work. Acknowledgment: We 
would like to thank J.L. Cardy for useful comments. 
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